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Abstract.  Aircraft  reconnaissance  data  collected  during  the 
Tropical  Cyclone  Structure  2008  field  campaign  are  used  to 
examine  further  kinematical,  dynamical,  and  thermodynam¬ 
ical  aspects  of  the  genesis  of  Typhoon  Nuri.  Data  from  the 
first  two  missions  into  the  pre-Nuri  disturbance  document 
the  transition  from  a  tropical  wave  to  a  tropical  depression. 
Dropwindsonde-derived  tangential  wind  profiles  at  several 
radii  from  the  low-level  circulation  center  indicate  that  the 
magnitude  of  low-level  circulation  increases  and  that  the  cor¬ 
responding  tangential  velocity  maximum  moves  inward  from 
the  first  to  second  reconnaissance  mission.  To  compliment 
these  findings,  a  three-dimensional  variational  analysis  in¬ 
corporating  both  dropwindsonde  and  aircraft  Doppler  radar 
data  is  conducted.  These  data  are  used  to  perform  circula¬ 
tion  tendency  calculations  at  multiple  distances  from  the  low- 
level  circulation  center.  The  results  demonstrate  a  net  spin- 
up  of  the  system-scale  circulation  in  the  low  levels  near  the 
center  and  in  the  outer  regions  of  the  recirculating  Kelvin 
cat’s  eye  circulation.  In  these  regions,  the  spin-up  tendency 
due  to  the  influx  of  cyclonic  absolute  vorticity  exceeds  the 
frictional  spin-down  tendency  for  both  Nuri  missions.  The 
system-scale  spin-up  is  found  to  be  accompanied  by  areas  of 
low-level  vorticity  concentration  through  vortex-tube  stretch¬ 
ing  associated  with  cumulus  convection.  The  areal  coverage 
and  intensity  of  these  high-vorticity  regions  increase  between 
the  first  and  second  Nuri  missions. 

The  findings  of  this  study  are  consistent  in  some  re¬ 
spects  to  the  Nuri  observational  analysis  carried  out  by  Ray¬ 
mond  and  Lopez-Carrillo  (201 1),  but  differ  in  their  suggested 
key  results  and  related  scientific  implications  that  the  pre- 


Nuri  disturbance  was  spinning  down  in  the  planetary  bound¬ 
ary  layer  on  the  first  day  of  observations.  The  findings  herein 
strongly  support  a  recent  tropical  cyclogenesis  model  posit¬ 
ing  that  the  Kelvin  cat’s  eye  circulation  of  the  parent  wave¬ 
like  disturbance  provides  a  favorable  environment  for  con¬ 
vective  vorticity  organization  and  low-level  spin-up  on  the 
mesoscale. 


1  Introduction 

Although  advances  in  the  understanding  and  prediction  of 
some  aspects  of  mature  tropical  cyclones  have  improved  over 
recent  years  (AMS  Statement,  2007),  there  remain  signif¬ 
icant  challenges  to  fully  understand  the  problem  of  tropi¬ 
cal  cyclogenesis.  Two  significant  hurdles  that  must  be  over¬ 
come  are  (i)  gaining  a  complete  understanding  of  the  com¬ 
plex  multi-scale  interactions  that  occur  during  the  transition 
of  a  synoptic-scale  precursor  to  a  mesoscale  vortex  and  (ii) 
addressing  the  lack  of  in  situ  observations  in  areas  where 
seedling  disturbances  for  tropical  cyclones  originate  and  in¬ 
tensify. 

In  their  observational  study,  Dunkerton  et  al.  (2009,  here¬ 
after  DMW09)  developed  a  new  tropical  cyclogenesis  model 
that  outlines  the  dynamics  and  thermodynamics  of  the  gen¬ 
esis  sequence  on  the  synoptic,  meso-,  and  convective  scales. 
This  model  was  presented  for  disturbances  originating  from 
African  easterly  waves  or  other  tropical  waves  and  related 
hydrodynamic  instabilities  in  a  mean  shear  flow.  Using  three 
independent  data  sets,  European  Centre  for  Medium- Range 
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Weather  Forecasts  Reanalysis  data  (ECMWF  ERA-40), 
Tropical  Rainfall  Measurement  Mission  (TRMM)  3B42  3- 
hourly  precipitation  and  best- track  data  from  the  National 
Flurricane  Center  (NFIC),  the  Kelvin  cat’s  eye  within  the  crit¬ 
ical  layer  of  a  tropical  easterly  wave  was  shown  to  be  the  pre¬ 
ferred  region  of  storm  formation,  more  precisely  located  near 
the  intersection  of  the  wave  critical  latitude  and  trough  axis. 
The  wave  critical  layer  is  the  region  surrounding  a  nonlin¬ 
ear  wave  critical  latitude  or  surface  in  latitudinal  shear  flow 
where  the  wave’s  phase  speed  equals  the  mean  flow  velocity. 
In  the  enclosed  Kelvin  cat’s  eye,  particles  are  trapped  and 
recirculate,  rather  than  being  swept  one  way  or  the  other  by 
the  surrounding  latitudinal  shear.  DMW09  articulated,  in  the 
form  of  scientific  hypotheses,  three  main  ways  that  the  recir¬ 
culating  cat’s  eye  region  of  a  tropical  wave  provide  a  favor¬ 
able  environment  for  tropical  cyclogenesis:  (HI)  wave  break¬ 
ing  or  roll-up  of  the  cyclonic  vorticity  and  moisture  near  the 
critical  surface  in  the  lower  troposphere  provides  the  moist 
vorticity  seedlings  and  a  favorable  environment  for  vorticity 
aggregation  leading  to  tropical  cyclone  formation;  (H2)  the 
cat’s  eye  is  a  region  of  quasi-closed  Lagrangian  circulation, 
and  air  is  repeatedly  moistened  by  convection  and  protected 
to  some  degree  from  dry  air  intrusion,  which  favors  a  pre¬ 
dominantly  convective  type  of  heating  profile;  (H3)  the  par¬ 
ent  wave  is  maintained  and  possibly  enhanced  by  diabatically 
amplified  mesoscale  vortices  within  the  cat’s  eye.  The  com¬ 
bination  of  the  foregoing  hypotheses  H1-H3  was  labeled  the 
marsupial  paradigm  by  DMW09  insofar  as  the  cyclogenesis 
sequence  is  likened  to  the  development  of  a  marsupial  infant 
in  its  mother’s  pouch.  The  “juvenile”  proto-vortex  is  carried 
along  and  protected  by  the  “mother”  wave  until  it  is  strength¬ 
ened  into  a  self-sustaining  entity.  The  cat’s  eye  within  the 
wave  critical  layer  is  thus  dubbed  the  “wave’s  pouch”  or  sim¬ 
ply,  “pouch”. 

Data  collected  during  two  recent  field  campaigns  have  pro¬ 
vided  in  situ  observations  to  test  the  foregoing  hypotheses 
under  various  atmospheric  conditions  and  in  different  tropi¬ 
cal  cyclone  basins.  The  THORPEX  (The  Observing  System 
Research  and  Predictability  Experiment)  Pacific  Regional 
Campaign  (TPARC)  and  Tropical  Cyclone  Stmcture  2008 
(TCS-08)  experiment  took  place  in  the  western  North  Pa¬ 
cific  basin  in  August  and  September  of  2008  (Elsberry  and 
Harr,  2008).  Data  collected  from  these  campaigns  provided 
the  first  opportunity  to  test  the  applicability  of  the  new  trop¬ 
ical  cyclogenesis  model  in  the  western  North  Pacific  basin 
for  the  case  of  westward  propagating  disturbances.  The  Pre- 
Depression  Investigation  of  Cloud  systems  in  the  Tropics 
(PREDICT)  campaign  took  place  in  the  Atlantic  basin  in 
2010  (Montgomery  et  al.,  2012).  This  campaign  focused  ex¬ 
clusively  on  collecting  data  in  the  pre-formation  stage  of  po¬ 
tential  tropical  cyclones.  The  tropical  cyclogenesis  cases  of 
Karl,  Matthew,  and  Gaston  have  been  examined  in  the  con¬ 
text  of  this  new  tropical  cyclone  model  (Smith  and  Mont¬ 
gomery,  2012;  Rutherford  and  Montgomery,  2012;  Wang, 


2012).  Analyses  of  data  from  these  two  campaigns  have  pro¬ 
vided  new  insight  into  the  tropical  cyclogenesis  problem. 

The  case  of  Typhoon  Nuri  was  one  of  the  more  robustly 
sampled  disturbances  among  these  campaigns  as  aircraft  re¬ 
connaissance  missions  flew  through  the  disturbance  on  four 
consecutive  days,  allowing  observation  of  the  storm  during 
the  tropical  wave,  tropical  depression,  tropical  storm,  and  ty¬ 
phoon  stages.  Montgomery  et  al.  (2010a)  used  observations 
and  global  model  analyses  to  examine  the  formation  of  Ty¬ 
phoon  Nuri  on  the  synoptic  and  meso-a  scales.  They  showed 
that  the  synoptic-scale  precursor  disturbance  was  a  tropical 
easterly  wave  with  maximum  amplitude  in  the  lower  tropo¬ 
sphere  (850  hPa  and  below).  The  characteristics  of  the  wave 
including  propagation  speed  and  direction,  low-level  vortic¬ 
ity  maximum,  and  distribution  of  precipitation  were  shown 
to  be  consistent  with  previous  observational  studies  of  east¬ 
erly  waves  in  the  western  North  Pacific  basin  dating  back  to 
the  1970s  (e.g.,  Chang,  1970;  Reed  and  Recker,  1971).  The 
wave-like  characteristics  of  this  disturbance  provided  the  op¬ 
portunity  for  Montgomery  et  al.  (2010a)  to  test  the  applica¬ 
bility  of  the  marsupial  paradigm  in  the  western  North  Pa¬ 
cific  basin.  Montgomery  et  al.  (2010a)  showed  that  a  trop¬ 
ical  wave  in  the  western  North  Pacific  basin  has  the  ability 
to  form  a  protective  wave  pouch  that  is  more  moist  than  the 
surrounding  environment.  Additionally,  they  showed  that  Ty¬ 
phoon  Nuri’s  wave  pouch  was  favorable  for  convection  and 
that  the  “sweet  spot”  proved  to  be  the  favored  region  for  low- 
level  vorticity  aggregation  and  growth  of  the  pro  to- vortex. 

While  Montgomery  et  al.  (2010a)  provided  an  overview  of 
the  tropical  cyclogenesis  sequence  on  larger  scales,  they  did 
not  address  mesoscale  organization  within  the  wave  pouch. 
A  complimentary  study  by  Raymond  and  Lopez-Carrillo 
(2011,  hereafter  RLll)  used  dropwindsonde  observations 
and  airborne  Doppler  radar  retrievals  to  examine  the  circula¬ 
tion  budget  of  Typhoon  Nuri.  The  analyses  of  RLl  1  showed 
that  for  the  first  research  flight  (hereafter  Nuri  1),  there  is 
an  absolute  circulation  maximum  from  the  surface  to  4  km 
that  is  nearly  constant  in  magnitude  with  height.  Contrarily, 
RLl  1  show  that  the  vertical  profile  of  absolute  circulation  in 
the  second  Nuri  flight  (hereafter  Nuri  2)  revealed  a  distinct 
maximum  at  the  5  km  level.  By  calculating  the  circulation 
tendency  over  an  area  roughly  enclosed  by  the  dropwind¬ 
sonde  locations  for  each  mission,  RLl  1  concluded  that  there 
was  a  net  spin-down  of  the  vortex  in  the  planetary  bound¬ 
ary  layer  (PBL)  during  Nuri  1  on  account  of  surface  friction 
dominating  the  spin-up  tendency  associated  with  the  influx 
of  cyclonic  absolute  vorticity  into  the  analysis  region. 

RLll  and  Raymond  et  al.  (2011)  reasoned  that  since 
Nuri  1  was  spinning  down  in  the  lower  levels,  a  thermody¬ 
namic  mechanism  was  required  to  reinitiate  low-level  spin- 
up.  The  basis  for  this  hypothesis  can  be  traced  back  to  a 
set  of  cumulus  ensemble  numerical  experiments  in  a  two- 
dimensional,  non-rotating  equatorial-like  environment  per¬ 
formed  by  Raymond  and  Sessions  (2007).  In  that  study,  the 
authors  found  that  increased  stabilization  of  the  atmosphere 
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(through  the  generation  of  a  mid-level  warm  temperature  per¬ 
turbation  coupled  with  the  generation  of  a  cold  temperature 
perturbation  in  the  lower  troposphere)  resulted  in  a  concen¬ 
trated  inflow  in  a  shallow  layer  in  the  lower  troposphere. 
They  found  also  that  the  thermal  stabilization  was  more  ef¬ 
fective  than  the  moistening  of  a  column  in  fostering  low-level 
convergence.  The  implication  is  that  if  realistic  values  of  am¬ 
bient  rotation  associated  with  a  tropical  wave  or  monsoon 
trough  were  included,  this  convergence  would  import  suffi¬ 
cient  absolute  vorticity  to  overcome  that  lost  to  the  surface 
by  friction.  This  thermodynamic  configuration  was  argued  to 
be  an  essential  ingredient  in  the  spin-up  of  a  tropical  cyclone 
by  altering  the  vertical  mass  flux^  profile  and  suppressing 
the  lateral  export  of  moist  entropy  (Raymond  et  al.,  2011).  In 
the  case  of  Nuri  1 ,  RLl  1  showed  that  the  area-averaged  ver¬ 
tical  mass  flux  maximum  was  located  near  10  km,  while  the 
lower  levels  (below  4  km)  exhibited  a  near-zero  vertical  mass 
flux  profile  with  a  local  region  of  negative  values,  presum¬ 
ably  in  association  with  mesoscale  convective  downdrafts. 
In  Nuri  2,  the  vertical  mass  flux  maximum  was  found  near 
5  km.  They  postulated  that  as  the  height  of  the  vertical  mass 
flux  maximum  is  lowered,  the  low-level  convergence  would 
be  enhanced  and  thus,  the  processes  described  above  would 
increase  the  low-level  vorticity. 

While  RLll  and  Raymond  et  al.  (2011)  offer  a  useful 
and  insightful  compliment  to  the  large-scale  overview  pre¬ 
sented  in  Montgomery  et  al.  (2010a),  there  remain  funda¬ 
mental  questions  regarding  the  mesoscale  organization  of 
Typhoon  Nuri  that  need  to  be  addressed.  Montgomery  and 
Smith  (2012)  first  questioned  the  conclusions  presented  in 
RLll  and  Raymond  et  al.  (2011).  Montgomery  and  Smith 
(2011)  used  dropwindsonde  data  from  the  Nuri  missions  to 
examine  the  thermodynamic  characteristics  of  the  recircu¬ 
lating  pouch  region  of  the  wave-like  disturbance  in  which 
Nuri  formed.  They  showed  that,  on  the  scale  of  the  pouch 
region,  the  cooling  in  the  lower  troposphere  between  Nuri  1 
and  2  was  negligible  and  the  warming  in  the  upper  tropo¬ 
sphere  was  less  than  1  K,  raising  significant  questions  con¬ 
cerning  the  necessity  of  the  thermodynamic  control  postu¬ 
lated  by  Raymond  et  al.  (2011).  Along  these  lines,  Wang 
(2012)  suggested  that  the  evolution  of  circulation  dynam¬ 
ics  and  thermal  stabilization  differ  on  different  spatial  scales 
within  the  wave  pouch.  In  the  inner  region,  where  convec¬ 
tive  processes  dominate,  Wang  (2012)  showed  that  low-level 
spin-up  occurs  and  thermal  stabilization  is  not  present.  Ad¬ 
ditionally,  Wang  (2012)  showed  that  the  outer  pouch  region 
is  characterized  by  weak  thermal  stabilization  and  is  accom¬ 
panied  by  low-level  spin-down. 


*  We  note  that  the  vertical  mass  flux  presented  in  RLl  1  is  actu¬ 
ally  an  area  integrated  quantity  (with  units  of  kgs“^)  as  opposed 
to  a  true  flux  quantity  (with  units  of  kg  s“  ^  m“^) .  In  this  paper,  we 
calculate  the  area  integrated  vertical  mass  flux  as  in  RLll  and  re¬ 
fer  to  the  quantity  as  vertical  mass  flux  (as  in  RLll)  to  maintain 
consistency. 


Additionally,  Smith  and  Montgomery  (2012)  and  Mont¬ 
gomery  and  Smith  (2012)  argued  that  the  tropical  cyclogene¬ 
sis  model  derived  from  the  non-rotating,  radiative  convective 
equilibrium  simulations  presented  in  Raymond  and  Sessions 
(2007)  is  not  strictly  appropriate  for  off-equatorial  tropical 
disturbances,  such  as  the  case  of  the  pre-Typhoon  Nuri  dis¬ 
turbance.  Finally,  differences  in  data  presented  in  RLll  and 
Raymond  et  al.  (2011)  lead  to  potential  confusion  on  the  dy¬ 
namical  mechanism  of  low-level  spin-up.  Specifically,  Ray¬ 
mond  et  al.  (2011)  showed  that  the  area- averaged  low-level 
absolute  vorticity  increased  from  Nuri  1  to  2,  contrary  to  their 
previous  finding  in  RLl  1  that  there  was  low-level  spin-down. 
The  area  in  which  the  average  low-level  absolute  vorticity 
was  calculated  differs  slightly  between  RLll  and  Raymond 
et  al.  (2011),  suggesting  that  part  of  the  discrepancy  may  be 
due  to  the  different  averaging  areas.  Similarly,  Raymond  et 
al.  (2011,  their  Fig.  9)  showed  that  the  average  vorticity  ten¬ 
dency  was  positive  throughout  the  lower  troposphere  except 
at  the  surface  (where  it  is  —0.001  ks“^  dy~^  or  nearly  zero). 
These  results  again  differ  from  RLll^.  Given  these  appar¬ 
ent  inconsistencies  and  the  importance  of  the  scientific  is¬ 
sues  summarized  above,  we  believe  there  is  a  need  for  an 
independent  assessment  of  the  dynamics  of  Nuri’s  genesis. 
The  objective  of  this  study,  then,  is  to  answer  the  following 
questions: 

(i)  Was  there  spin-down  on  the  system  scale  in  the  lower 
levels  of  the  developing  Typhoon  Nuri  disturbance  be¬ 
tween  Nuri  1  and  2  as  suggested  by  RLl  1  ? 

(ii)  Is  the  thermodynamic  control  hypothesis  put  forth  by 
Raymond  et  al.  (2011)  a  necessary  condition  for  tropical 
cyclogenesis  in  the  case  of  Typhoon  Nuri? 

(iii)  How  does  the  chosen  area  of  integration  for  circulation 
budget  calculations  affect  the  results  and  the  conclu¬ 
sions  drawn  from  these  analyses? 

The  outline  of  this  paper  is  as  follows:  Sect.  2  discusses 
the  data  and  methodology  used  in  our  analysis.  Section  3 
presents  the  results  of  our  analysis  for  Nuri  1  and  2.  Section  4 
provides  a  discussion  of  our  results  and  their  implications  on 
our  understanding  of  Nuri’s  formation.  Section  5  presents  the 
conclusions. 

2  Data  and  methodology 

2.1  Research  flights  and  data  collected 

The  precursor  easterly  wave  that  contributed  to  the  formation 
of  Typhoon  Nuri  was  first  identified  in  operational  forecast 

^The  differences  in  wave  propagation  speed  used  in  RLl  1 
(— 5  ms“^)  versus  Raymond  et  al.  (2011)  (—7  ms~^)  may  have  led 
to  the  small  differences  in  magnitude  of  circulation  (vorticity)  ten¬ 
dency. 
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Nuri  1  Nuri  2 


Figure  1.  Electra-Doppler  Radar  reflectivity  composite  (in  dBZ)  at  3  km  altitude  for  Nuri  1  (left)  and  Nuri  2  (right)  superimposed  with 
dropwindsonde  locations.  The  dropwindsondes  from  the  USAF  C130  (filled  circles)  and  NRL  P3  (filled  squares)  have  been  translated  to 
their  00:00  UTC  positions  using  a  (westward)  phase  speed  of  — 7  m  s~  * .  The  red  circles  and  boxes  indicate  dropwindsondes  with  no  retrieved 
wind  data.  The  square  boxes  are  at  one-degree  increments,  increasing  from  one  degree  to  four  degrees,  and  represent  the  domain  used  in 
several  of  the  forthcoming  analyses.  The  boxes  are  centered  on  the  sweet  spot  position  at  1 .5  km  altitude  within  the  pre-Nuri  wave  as  derived 
from  the  25  km  SAMURAI  analyses.  The  abscissa  is  longitude  and  the  ordinate  is  latitude. 


models  on  the  morning  of  12  August  2008^.  The  first  TCS- 
08  research  flight  was  a  joint  United  States  Air  Force  (USAF) 
Cl 30  and  Naval  Research  Laboratory  (NRL)  P3  flight  into 
the  developing  wave  on  15-16  August  2008.  The  second  re¬ 
search  flight  was  flown  on  16-17  August  2008  and  was  also 
a  joint  USAF  C130,  NRL  P3  mission.  The  Joint  Typhoon 
Warning  Center  (JTWC)  upgraded  the  disturbance  to  a  tropi¬ 
cal  depression  (TD)  at  12:00 UTC,  16  August  2008, between 
the  first  and  second  research  flights.  The  third  research  flight 
included  both  aircraft  and  occurred  between  17-18  August, 
after  JTWC  upgraded  the  disturbance  to  a  tropical  storm.  The 
final  research  mission  was  a  USAF  C130-only  mission  on 
18-19  August,  when  the  disturbance  was  a  typhoon.  Follow¬ 
ing  RLl  1  and  consistent  with  the  terminology  defined  in  the 
Introduction,  we  refer  to  the  state  of  Nuri  observed  during 
these  four  observation  times  as  Nuri  1, 2,  3  and  4,  respec¬ 
tively.  Since  this  paper  is  focused  on  the  tropical  cyclogene¬ 
sis  sequence  and  not  on  the  subsequent  intensification  of  the 
tropical  cyclone,  we  will  focus  primarily  on  the  data  from  the 
first  two  missions^. 

The  two  types  of  data  we  use  in  this  analysis  are  drop¬ 
windsonde  observations  from  both  the  USAF  Cl 30  and  NRL 
P3  and  Electra-Doppler  Radar  (ELDORA)  data  from  the 
NRL  P3.  The  USAF  Cl 30  released  the  dropwindsondes  from 
heights  of  up  to  9.4  km  while  the  NRL  P3  typically  operated 


^Montgomery  et  al.  (2010a)  postulated  that  the  wave  originated 
on  7  August  in  the  central  Pacific. 

^The  exception  is  our  Fig.  9,  where  for  completeness  and  to 
compare  our  specific  findings  with  RLll,  we  present  the  relative 
frequency  distribution  of  vorticity  for  Nuri  1 , 2  and  3. 


at  lower  elevations  (between  2.6  and  3.4km).  The  NRL  P3 
flew  at  higher  elevations  during  the  end  of  the  first  Nuri  mis¬ 
sion  (approximately  7.3  km)  to  collect  more  extensive  drop¬ 
windsonde  profiles.  The  ELDORA  is  an  airborne  dual-beam 
meteorological  radar  with  an  unambiguous  range  of  75  km. 
Doppler  radar  data  were  collected  in  150  km  swaths  centered 
on  the  P3  track  (RLl  1).  The  radar  data  were  edited  using  the 
quality  control  techniques  described  in  Bell  et  al.  (2013).  Fig¬ 
ure  1  illustrates  the  spatial  coverage  of  the  Cl 30  dropwind¬ 
sonde  observations  (filled  circles),  NRL  P3  dropwindsondes 
(filled  squares),  and  ELDORA  coverage  (as  indicted  by  the 
3  km  altitude  reflectivity  return,  dBZ)  for  Nuri  1  (left)  and 
Nuri  2  (right).  In  all  of  the  3-D  variational  (3-D- VAR)  anal¬ 
yses  presented,  these  data  are  translated  to  their  00:00  UTC 
positions  based  on  a  zonal  phase  speed  of  — 7  m  s“' .  The  thin 
boxes  in  Fig.  1  are  centered  on  the  sweet  spot  position  at 
1 .5  km  altitude  out  to  four  degrees  length  in  one-degree  in¬ 
crements.  We  choose  the  1 .5  km  altitude  sweet  spot  position 
for  these  calculations  because  the  maximum  amplitude  of  the 
pre-Nuri  wave  was  in  the  lower  troposphere  (Montgomery 
et  al.,  2010a)^.  The  boxes  are  illustrative  in  nature  and  pro¬ 
vide  the  reader  with  a  sense  of  the  data  coverage  at  varying 
distances  from  the  sweet  spot.  The  forthcoming  circulation 
tendency  calculations  are  conducted  on  similar  boxes  at  0.2 
degree  increments  from  the  sweet  spot  center.  It  is  noted  that 


^The  choice  of  1.5  km  altitude  is  consistent  with  the  observed 
properties  of  easterly  waves  in  the  western  North  Pacific  basin  dat¬ 
ing  back  to  the  1970s,  which  found  that  the  maximum  amplitude  of 
easterly  waves  was  in  the  lower  troposphere,  at  or  below  850  hPa 
(Reed  and  Recker,  1971;  Chang,  1970). 


Atmos.  Chem.  Phys.,  14, 8795-8812, 2014 


www.atmos-chem-phys.net/14/8795/2014/ 


L.  L.  Lussier  III  et  al.:  The  genesis  of  Typhoon  Nuri  as  observed  during  the  TCS-08  field  experiment 


8799 


several  of  the  dropwindsondes  (indicated  by  red  coloring)  do 
not  contain  velocity  data  due  to  instrument  errors. 

To  provide  a  foundation  for  assessing  the  changes  in  low- 
level  spin-up  between  Nuri  1  and  2,  we  first  carry  out  a  cir¬ 
culation  analysis  using  only  dropwindsonde  data.  The  coarse 
resolution  and  irregular  spacing  of  the  dropwindsonde  data 
make  identifying  a  sweet  spot  position  challenging.  To  iden¬ 
tify  the  cyclonic  circulation  center  at  each  level,  we  fol¬ 
low  the  methodology  of  Davis  and  Ahijevych  (2012).  This 
method  identifies  the  maximum-weighted  average  tangen¬ 
tial  wind  within  a  three-degree  radius  circle  at  each  location 
within  the  dropwindsonde  domain.  We  perform  this  calcula¬ 
tion  at  5  hPa  intervals  in  the  vertical  and  the  horizontal  lo¬ 
cation  where  this  average  tangential  wind  is  maximized  is 
designated  as  the  cyclonic  circulation  center.  The  circulation 
centers  at  each  level  are  used  to  examine  the  vertical  tilt  of 
the  system.  However,  the  850  hPa  circulation  center  is  cho¬ 
sen  as  the  center  position  to  calculate  tangential  wind  at  all 
levels  for  Nuri  1  and  2. 

To  compliment  the  dropwindsonde-only  analysis,  EL- 
DORA  and  dropwindsonde  data  are  analyzed  together  in 
a  newly  developed  3-D-VAR  analysis  scheme  (Bell  et 
al.,  2012)  referred  to  as  SAMURAI  (Spline  Analysis  at 
Mesoscale  Using  Radar  and  Aircraft  Instrumentation).  The 
SAMURAI  analysis  uses  a  Bayesian  approach  to  determine 
the  maximum  likelihood  estimate  of  the  thermodynamic  and 
kinematic  data  from  the  observations  (Bell  et  al.,  2012). 
SAMURAI  is  similar  to  other  variational  schemes  that  min¬ 
imize  the  difference  between  observations  and  the  resulting 
analysis,  but  unlike  numerical  weather  prediction  data  assim¬ 
ilation  it  enforces  no  balance  or  physical  constraints  other 
than  mass  continuity.  The  analysis  scheme  uses  a  series  of  fi¬ 
nite  elements  in  the  form  of  cubic  B-splines  to  represent  the 
3-D  atmospheric  state.  The  spline  analysis  has  the  advantage 
of  near-analytic  derivatives  for  vorticity  and  divergence,  and 
precise  interpolation  to  the  observation  space  for  irregularly 
spaced  data  such  as  dropwindsonde  and  airborne  radar  data. 
SAMURAI  bears  some  similarity  to  the  3-D-VAR  software 
used  in  RLl  1  (Lopez-Carrillo  and  Raymond,  2011),  but  has 
a  distinctly  different  numerical  implementation  of  the  vari¬ 
ational  cost  function.  No  fundamental  differences  between 
the  wind  fields  presented  here  and  those  presented  in  Rll 
were  apparent,  suggesting  that  both  analysis  techniques  cap¬ 
ture  the  salient  features  of  Nuri,  but  differ  in  the  details. 

Several  SAMURAI  analyses  were  created  for  this  study. 
For  each  research  flight,  a  synoptic-scale  analysis,  with  hori¬ 
zontal  grid  spacing  of  25  km,  and  a  mesoscale  analysis,  with 
horizontal  grid  spacing  of  10km,  were  created.  A  Gaussian 
low-pass  filter  was  applied  to  the  analyses  to  spread  infor¬ 
mation  from  the  observations  and  reduce  noise  in  the  wind 
fields.  The  horizontal  filter  lengths  were  125  and  50  km  for 
the  synoptic-scale  and  mesoscale  analyses,  respectively.  All 
analyses  have  0.5  km  vertical  grid  spacing  from  the  surface 
to  an  upper  height  of  16.5  km,  with  a  1km  vertical  filter 
length.  In  the  forthcoming  analyses,  we  only  present  data  up 


to  12  km  in  the  vertical  due  to  the  high-elevation  sampling 
limitations  of  the  ELDORA.  The  main  purpose  of  the  25  km 
analyses  is  to  identify  the  sweet  spot  position  during  Nuri  1 
and  2  as  it  is  hypothesized  to  be  the  region  for  upscale  vor¬ 
ticity  growth  and  organization  (DMW09,  Montgomery  et  al., 
2010b).  The  10km  analyses  are  used  primarily  to  construct 
the  kinematic  and  dynamic  variables  for  the  vorticity  analy¬ 
sis. 

2.2  Vorticity  and  circulation  dynamics 

The  vorticity  dynamics  are  diagnosed  using  the  10  km 
SAMURAI  analyses  following  RLl  1 .  Following  Haynes  and 
McIntyre  (1987)  and  RLll,  the  vertical  vorticity  equation 
may  be  rewritten  into  the  flux  form  as  follows: 

^  =  -Vh-z  +  k-Vh0xVhn,  (1) 

at 

where  t]  is  the  absolute  vertical  vorticity,  the  subscript  h 
represents  the  horizontal  vector  quantity,  and  6  is  poten¬ 
tial  temperature.  The  Exner  function  is  defined  as  follows: 
n  =  (^)^‘*/'^p,  where  p  is  the  pressure,  po  is  the  pressure  at 
a  reference  level  (typically  1000  hPa),  is  the  gas  constant 
for  dry  air,  and  Cp  is  the  heat  capacity  of  dry  air  at  constant 
pressure.  The  vector  vorticity  flux,  Z,  is  defined  as  follows: 

Z  =  Zi+Z2  +  Z3  =  ri^Vh-VhVz  +  k><  F,  (2) 

where  the  z-subscript  represents  the  vertical  component  of 
the  vector  and  the  vector  F  represents  horizontal  frictional 
and  other  sub-grid-scale  forces.  If  we  assume  that  F  can  be 
approximated  by  friction  only,  it  can  be  parameterized  in  a 
simple  and  plausible  way  (following  RLll)  as  follows: 

F  =-exp{-z/zs)/zs,  (3) 

P 

where  p  is  air  density,  z  is  height,  Zs  is  a  characteristic  bound¬ 
ary  layer  height  ( 1 .25  km)  and  the  surface  stress  t  is  defined 
as  follows: 

T  = -pblCd  IUblI  (4) 

where  the  subscript  BE  represents  conditions  in  the  boundary 
layer,  U  is  the  horizontal  wind  and  Cd  is  the  surface  drag 
coefficient,  defined  as  follows: 

Cd  =  (1 +0.028  IUblDx  10-3  (5^ 

Now,  integrating  Eq.  (1)  over  a  horizontal  surface,  and  ne¬ 
glecting  the  small  contribution  from  the  baroclinic  vorticity 
generation  term,  results  in  the  circulation  tendency  equation: 

JEa  f 

— ^  =  /  W  -ZdA,  (6) 

dt  J 

A 
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where  Fa  is  the  absolute  circulation  and  A  is  the  area  over 
which  the  integration  is  performed  (in  this  work,  a  square 
box).  By  applying  Gauss’s  divergence  theorem,  the  circula¬ 
tion  tendency  can  also  be  calculated  using  the  line  integral 
form: 


dt 


v„r]^dl  + 


rinVzdl 


Ftdl, 


(7) 


where  dl  is  the  positive  line  element  along  the  perimeter 
of  the  integration  area,  the  subscript  n  denotes  the  outward 
normal  component  of  the  variable  along  the  circuit,  and  the 
subscript  t  denotes  the  tangential  component  of  the  vari¬ 
able  in  the  sense  of  the  circuit.  The  integration  is  taken  in 
a  counterclockwise  sense  in  accordance  with  the  right-hand 
mle.  Equation  (7)  shows  that  changes  to  the  absolute  cir¬ 
culation  occur  through  (i)  convergence  of  absolute  vorticity 
(first  term);  (ii)  vortex  tilting-like  term  (second  term);  and 
iii)  the  frictional  spin-down  tendency  (third  term).  RLll  ar¬ 
gue  that  the  spin-up  of  a  vortex  in  the  lower  levels  can  oc¬ 
cur  only  if  the  spin-up  from  convergence  of  absolute  vor¬ 
ticity  exceeds  the  frictional  spin -down  tendency®.  Davis  and 
Galameau  (2011)  show  that  the  convergence  term  can  be  bro¬ 
ken  down  into  “mean”  and  “eddy”  contributions  along  the 
circuit: 


—  ®  v„rizdl  =  —rjSA  —  ®  rj'v'  ■  hdl, 


(8) 


where  S  is  the  horizontal  divergence,  the  overbars  represent 
the  mean  over  the  perimeter  of  the  circuit,  the  primes  indi¬ 
cate  perturbations  from  this  mean,  the  tildes  indicate  an  areal 
average  over  the  circuit,  and  v  is  shorthand  for  the  horizontal 
velocity  vector.  The  first  term  on  the  RHS  of  Eq.  (8)  repre¬ 
sents  vortex-tube  stretching  averaged  over  the  area  contained 
within  the  circuit,  while  the  second  term  in  this  equation  rep¬ 
resents  the  change  in  the  absolute  circulation  owed  to  hor¬ 
izontal  eddy  fluxes  of  vertical  vorticity  into  and  out  of  the 
integration  area. 

^Several  previous  studies  (Willoughby,  1995;  Schubert  and  Al- 
worth,  1987;  Marin  et  al.,  2009)  share  a  similar  view  with  RLl  1  and 
Raymond  et  al.  (201 1)  in  that  they  suggest  that  friction  within  the 
boundary  layer  acts  solely  as  a  detriment  to  spin-up  on  the  system- 
scale  circulation.  We  do  not  dispute  the  strict  statement  that  spin- 
up  in  the  boundary  layer  only  occurs  when  convergence  of  abso¬ 
lute  vorticity  is  larger  than  the  tendency  to  spin  down  the  flow  due 
to  friction.  It  is  not  true,  however,  that  the  dynamics  of  the  fric¬ 
tional  boundary  layer  act  solely  as  a  detriment  to  spin-up.  Smith  et 
al.  (2009),  Sanger  et  al.  (2013),  and  Montgomery  and  Smith  (2014) 
illustrate  that  within  the  inner-core  region  of  a  developing  tropical 
cyclone  vortex,  the  frictional  inflow  acquires  relatively  significant 
values  compared  to  the  tangential  component.  These  latter  studies 
demonstrate  a  dynamical  role  of  the  boundary  layer,  wherein  air 
parcels  may  reach  small  radii  quickly  (minimizing  the  loss  of  abso¬ 
lute  angular  momentum  (defined  in  Sect.  3)  during  spiral  circuits) 
and  therefore  acquiring  large  v  that  may  be  larger  than  the  v  in  the 
interior  flow  (and  the  gradient  value  thereof)  above  the  boundary 
layer. 


One  of  the  objectives  of  this  work  is  to  quantify  how  the 
circulation  dynamics  vary  within  the  wave  pouch  at  differ¬ 
ent  distances  from  the  sweet  spot.  With  this  objective  in 
mind,  the  SAMURAI  analyses  are  used  to  calculate  terms 
from  the  circulation  tendency  equation  (Eq.  7)  and  other  vari¬ 
able  properties  in  0.2  degree  length  boxes  centered  on  the 
1.5  km  altitude  sweet  spot  position.  The  advantages  of  this 
analysis  methodology  are  (i)  it  is  a  systematic  way  to  ana¬ 
lyze  the  kinematics  and  dynamics  at  varying  distances  from 
the  sweet  spot  position;  (ii)  it  allows  for  a  straightforward 
comparison  of  the  various  spin-up  contributions  between  re¬ 
search  flights;  (iii)  the  methodology  is  optimally  suited  to 
observe  and  quantify  the  development  that  ensues  around  the 
sweet  spot,  the  so-called  “attractor  point”  predicted  by  the 
new  cyclogenesis  model  (Montgomery  et  al.,  2012);  and  (iv) 
anchoring  our  analysis  to  the  sweet  spot  of  the  parent  distur¬ 
bance  provides  useful  information  on  the  vertical  alignment 
of  the  vortex  and  depth  of  the  wave  pouch.  The  second  and 
third  points  are  especially  important  as  there  is  considerable 
variability  in  the  areal  coverage  of  the  aircraft  data  during  the 
Nuri  missions  (Fig.  1). 

The  main  limitation  of  this  type  of  analysis  is  that  the  ob¬ 
served  data  are  not  distributed  uniformly  over  the  defined 
boxes.  This  results  in  data  being  unavailable  in  small  regions 
in  some  of  the  outer  boxes.  To  address  this  issue,  in  all  of 
the  calculations  performed  in  this  paper,  regions  outside  of 
the  observational  boundaries  are  set  to  missing  values.  Since 
these  issues  arise  either  at  the  outer  edge,  or  outside  of  the 
wave  pouch,  we  believe  that  the  impact  on  the  analysis  pre¬ 
sented  herein  is  minimal.  Non-uniform  data  distribution  is 
unfortunately  a  limitation  in  any  observational  study  and  we 
feel  that  the  advantages  gained  in  this  type  of  analysis  far 
outweigh  the  disadvantages. 


3  Results 

3.1  Dropwindsonde-only  analysis 

Dropwindsonde  analysis  from  the  first  research  flight  into 
the  pre-Nuri  (Nuri  1)  disturbance  is  presented  in  Fig.  2.  This 
analysis  depicts  wind  barbs  in  the  co-moving  frame  at  stan¬ 
dard  pressure  levels  up  to  500  hPa.  The  black  circles  indi¬ 
cate  the  position  of  the  850  hPa  circulation  center  and  the 
blue  triangles  indicate  the  position  of  the  circulation  center 
at  the  indicated  level.  The  low-level  (at  and  below  850  hPa) 
wind  fields  indicate  a  closed,  cyclonic  circulation  in  the  co¬ 
moving  frame  (Fig.  2,  top).  The  850  hPa  circulation  cen¬ 
ter  is  located  near  14.7°  N,  146.8°  E  and  is  nearly  vertically 
aligned  down  to  the  925  hPa  level.  At  700  and  500  hPa,  there 
are  cyclonic  circulation  centers  that  tilt  eastward  with  height. 
At  500  hPa,  the  cyclonic  circulation  center  is  located  near 
14.7°  N,  149.4°  E,  a  2.6  degree  eastward  tilt  with  increasing 
height. 
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Figure  2.  Four-panel  plot  of  wind  barbs  from  the  dropwindsonde 
observations  in  the  co-moving  frame  for  Nuri  1 .  Levels  plotted  are 
925  hPa  (a),  850  hPa  (b),  700  hPa  (c),  and  500  hPa  (d).  In  this  figure 
and  Fig.  3,  the  winds  are  in  standard  convention  for  wind  barbs 
(kt  -  knot),  rather  than  the  SI  units  used  in  the  remained  of  the 
paper.  A  full  wind  barb  denoted  10  kt  and  each  half  barb  is  5  kt.  The 
filled  circle  represents  the  850  hPa  dropwindsonde-derived  cyclonic 
circulation  center  and  the  filled  triangles  are  the  circulation  centers 
at  the  indicated  pressure  level.  The  abscissa  is  longitude  and  the 
ordinate  is  latitude. 


Figure  3.  Dropwindsonde  observations  from  the  Nuri  2  flight  (as  in 
Fig.  2). 

Figure  3  is  similar  to  Fig.  2,  but  for  Nuri  2.  The  850 hPa 
cyclonic  circulation  center  is  located  near  16.1°N,  140.3°E 
and  remains  nearly  vertically  aligned  down  to  925  hPa.  The 
circulation  center  has  become  more  longitudinally  aligned 


with  the  low-level  circulation  center  at  700  and  500  hPa, 
but  there  is  now  a  distinct  tilt  to  the  southeast  with  height. 
At  500  hPa,  the  cyclonic  circulation  center  is  located  near 
14.1°  N,  141 .1°  E,  a  2.2  degree  tilt  with  height.  Erom  a  qual¬ 
itative  point  of  view,  since  the  average  wind  speed  appears 
to  have  increased  at  all  levels  from  Nuri  1  to  2,  there  does 
not  appear  to  be  low-level  spin-down  between  these  two  mis¬ 
sions. 

Eigure  4  is  a  4-panel  plot  of  tangential  wind  profiles  de¬ 
rived  from  the  dropwindsonde  data  at  one-degree  radius  in¬ 
tervals  for  the  first  two  Nuri  missions.  Eor  each  dropwind¬ 
sonde,  the  tangential  wind  is  calculated  at  each  level  in  the 
vertical  and  then  binned  in  5  hPa  intervals.  An  average  tan¬ 
gential  wind  is  then  calculated  in  each  bin.  The  binning  is 
accomplished  to  smooth  high-resolution  variations  in  these 
data  and  as  an  attempt  to  limit  the  impact  of  missing  data 
when  direct  comparisons  are  made  between  dropwindsonde 
observations  at  particular  levels.  An  azimuthal  average  is 
then  performed  within  one-degree  radius  annuli  from  the 
850  hPa  dropwindsonde-derived  circulation  center. 

The  maximum  tangential  wind  for  Nuri  1  occurs  within 
the  2-3  degree  radius  annulus,  near  850 hPa.  The  inner  two 
annuli  indicate  areas  of  positive  tangential  wind  up  to  ap¬ 
proximately  600  hPa.  Above  this  level,  there  are  both  posi¬ 
tive  and  negative  values  of  tangential  wind,  likely  due  to  the 
tilt  of  the  vortex  with  height  discussed  above.  The  largest  in¬ 
crease  in  tangential  wind  between  Nuri  1  and  2  occurs  in  the 
1-2  degree  radius  annulus  at  600  hPa.  Eor  Nuri  2,  the  radius 
of  maximum  tangential  wind  is  less  well  defined,  as  all  radii 
beyond  one  degree  indicate  similar  wind  speeds.  The  1  de¬ 
gree  and  2-3  degree  annuli  show  slight  decreases  in  tangen¬ 
tial  wind  speed  from  Nuri  1  to  2  between  800-700  hPa.  The 
distribution  of  dropwindsonde  observations  for  Nuri  2  in  the 
outermost  radius  (Eig.  4d)  displays  a  strong  sampling  bias 
to  the  south  of  the  sweet  spot  position.  While  these  data  are 
shown  for  completeness,  it  is  difficult  to  draw  robust  con¬ 
clusions  with  respect  to  the  system-scale  circulation  in  the 
outer-most  radius  of  Nuri  2.  The  most  important  feature  evi¬ 
dent  in  these  profiles  is  that  there  is  an  increase  in  observed 
tangential  wind  speed  at  and  below  800  hPa  from  Nuri  1  to 
2  at  all  radii.  These  observations  suggest  that  spin-up  within 
and  just  above  the  boundary  layer  is  occurring  between  these 
two  missions. 

3.2  SAMURAI  analysis 

To  test  the  robustness  of  the  finding  that  low-level  spin-up 
is  occurring  between  the  first  two  missions,  the  SAMURAI 
analyses  are  examined  next.  Figure  5  shows  a  two-panel  plot 
of  the  Okubo- Weiss  (OW)  parameter,  streamlines  in  the  co¬ 
moving  frame,  and  isotachs  at  1 .5  and  5  km  altitude  from  the 
Nuri  1  mission  25  km  analysis.  Okubo- Weiss  is  defined  as 
OW  =  f  2  -  52  -  =  (V4  -  t/y)2  -  (Ux  -  Vy)2  -  (Ux  +  Uy)2 , 

where  f  is  relative  vorticity,  Si  and  S2  are  strain  deformation, 
{U,  V)  are  zonal  and  meridional  wind,  and  the  subscripts 


925  hPo 
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Figure  4.  Four-panel  plot  of  dropwindsonde-derived  azimuthally  averaged  tangential  wind  profiles  for  the  first  two  research  flights  into 
the  pre-Nuri  disturbance.  The  blue  curve  is  Nuri  1  and  the  red  curve  is  Nuri  2.  The  average  profiles  are  calculated  in  one-degree  radius 
annuli  from  the  dropwindsonde-derived  850  hPa  circulation  center  (filled  circles  in  Figs.  2  and  3)  for  each  research  flight.  The  insets  indicate 
dropwindsonde  positions  relative  to  the  low-level  circulation  center  for  each  annulus.  Each  circle  is  one  degree  radius.  The  blue  circles  are 
the  Nuri  1  USAF  C130  dropwindsonde  positions  and  the  blue  squares  are  the  Nuri  1  NRL  P3  dropwindsonde  positions.  The  circles  and 
squares  follow  the  same  convention  for  the  Nuri  2  dropwindsonde  positions. 


(x,  y)  indicate  partial  derivatives  in  the  zonal  and  merid¬ 
ional  directions,  respectively.  Significantly  positive  regions 
of  OW  indicate  areas  that  will  tend  to  be  immune  from  hori¬ 
zontal  shearing  deformation,  while  negative  regions  indicate 
areas  subject  to  shearing  deformation  (McWilliams,  1984; 
Weiss,  1991).  The  positions  of  the  wave  trough  (locus  of 
points  with  n  =  0)  and  critical  surface  (locus  of  points  with 
zero  relative  zonal  wind)  at  1.5km^  are  indicated  by  the 
thick  black  and  the  purple  lines,  respectively.  For  Nuri  1 ,  the 
1 .5  km  sweet  spot  position  is  located  near  14.6°  N,  146.3°  E. 
The  center  of  the  5  km  circulation  (as  determined  by  the  co¬ 
moving  streamlines)  is  near  12.9°  N,  147.9°  E.  This  indicates 
a  larger  tilt  to  the  southeast  with  height  than  observed  in 
the  dropwindsonde-only  data.  Although  the  exact  positions 
of  the  circulation  centers  are  slightly  different  from  those  in 
RLll,  the  overall  fiow  pattern  is  fairly  consistent  (cf.  RLll, 
p.  155).  At  1.5  km,  the  OW  maximum  is  centered  near  the 
sweet  spot  position  while  there  are  several  OW  maxima  at 

^While  the  dropwindsonde  data  are  displayed  on  pressure  levels 
in  the  vertical,  the  SAMURAI  data  are  displayed  in  height.  The 
1 .5  km  height  corresponds  to  approximately  the  850  hPa  level  and 
the  5  km  height  corresponds  to  approximately  the  550  hPa  level. 


5  km.  A  similar  analysis  for  Nuri  2  (Fig.  6)  indicates  a  1 .5  km 
sweet  spot  position  at  15.9° N,  140°  E.  The  5km  cyclonic 
circulation  has  now  become  longitudinally  aligned  with  the 
low-level  circulation  and  is  located  near  14.1°  N,  140.3°  E. 
The  OW  field  has  become  more  organized  and  concentrated 
near  the  cyclonic  circulation  center  at  both  levels. 

To  examine  the  observed  structure  of  vertical  vorticity  and 
horizontal  wind  within  the  wave  pouch  region.  Figs.  7  and  8 
show  two-panel  plots  of  absolute  vorticity  and  wind  vectors 
(co-moving)  from  the  10  km  SAMURAI  analyses  at  1 .5  and 
5  km  for  Nuri  1  and  2.  For  Nuri  1 ,  the  1 .5  km  mesoscale  cy¬ 
clonic  circulation  center  appears  north  of  the  25  km-analysis 
sweet  spot  position.  This  mesoscale  circulation  center  is  as¬ 
sociated  with  a  positive  vorticity  maximum  on  a  similar  hor¬ 
izontal  scale.  The  presence  of  mesoscale  circulation  centers 
that  are  not  co-located  with  the  sweet  spot  is  not  uncommon, 
as  multiple  studies  (Wang  et  al.,  2009;  Lussier  111,  2010) 
have  shown  the  presence  of  short-lived,  smaller-scale  circu¬ 
lations  within  the  larger  wave  pouch.  The  5  km  cyclonic  cir¬ 
culation  center  appears  to  be  southeast  of  the  low-level  cir¬ 
culation  center.  The  mesoscale  absolute  vorticity  maxima  at 
5  km  are  aligned  well  with  the  maxima  identified  in  the  25  km 
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Figure  5.  Two-panel  plot  of  Okubo-Weiss  (OW,  shaded)  and  streamlines  (co-moving)  for  Nuri  1  at  1 .5  km  (left)  and  5  km  (right)  height  from 
the  25  km  SAMURAI  analysis.  The  thick  purple  line  indicates  the  position  of  the  1 .5  km  critical  latitude  and  the  thick  black  line  indicates  the 
position  of  the  1.5  km  wave  trough.  The  dashed  yellow  lines  are  isotachs  contoured  at  2ms“*  intervals.  The  OW  is  only  shaded  in  regions 
with  cyclonic  relative  vorticity.  The  abscissa  is  longitude  and  the  ordinate  is  latitude. 
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Figure  6.  Two-panel  plot  of  Okubo-Weiss  (OW,  shaded)  and  streamlines  (co-moving)  for  Nuri  2  at  1 .5  km  (left)  and  5  km  (right)  height  from 
the  25  km  SAMURAI  analysis.  The  thick  purple  line  indicates  the  position  of  the  1 .5  km  critical  latitude  and  the  thick  black  line  indicates  the 
position  of  the  1.5  km  wave  trough.  The  dashed  yellow  lines  are  isotachs  contoured  at  2ms~^  intervals.  The  OW  is  only  shaded  in  regions 
with  cyclonic  relative  vorticity.  The  abscissa  is  longitude  and  the  ordinate  is  latitude. 


analysis.  The  overall  wind  field  and  vorticity  pattern  are  con¬ 
sistent  with  the  analysis  from  RLl  1 .  However,  the  magnitude 
of  the  absolute  vorticity  is  generally  weaker  in  our  analyses. 
The  Nuri  2  wind  field  exhibits  a  more  coherent  circulation 
than  observed  in  the  Nuri  1  analysis.  In  particular,  the  1 .5  km 
vorticity  maxima  have  become  more  consolidated  around  the 
sweet  spot  position.  There  are  localized  low-level  absolute 


vorticity  maxima  located  near  14.3°  N,  140.3°  E  and  16.8°  N, 
140.3°  E.  At  5  km  the  circulation  has  become  more  vertically 
aligned  with  height  in  the  longitudinal  direction  and  multiple 
vorticity  maxima  are  observed  within  the  wave  pouch. 

To  help  quantify  the  statistical  evolution  of  cyclonic  rel¬ 
ative  vorticity  that  takes  place  during  the  research  flights, 
Eig.  9  shows  a  4-panel  plot  of  the  relative  frequency  of 
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z=1.5  km  z  =  5  km 


Figure  7.  Two-panel  plot  of  absolute  vorticity  (shaded)  and  wind  vectors  (co-moving)  for  Nuri  1  at  1 .5  km  (left)  and  5  km  (right)  height 
from  the  10km  SAMURAI  analyses.  The  black  box  is  the  four-degree  length  box  centered  on  the  25  km  analysis  1.5  km  altitude  sweet  spot 
position  and  is  for  reference  only.  The  abscissa  is  longitude  and  the  ordinate  is  latitude. 
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Figure  8.  Two-panel  plot  of  absolute  vorticity  (shaded)  and  wind  vectors  (co-moving)  for  Nuri  2  at  1 .5  km  (left)  and  5  km  (right)  height 
from  the  10km  SAMURAI  analyses.  The  black  box  is  the  four-degree  length  box  centered  on  the  25  km  analysis  1.5  km  altitude  sweet  spot 
position  and  is  for  reference  only.  The  abscissa  is  longitude  and  the  ordinate  is  latitude. 


absolute  vorticity  for  the  first  three  Nuri  missions  derived 
from  the  10km  SAMURAI  analyses.  Frequency  distribu¬ 
tions  of  absolute  vorticity  within  0.5,  1,  1.5  and  2  degrees 
radius  of  the  1 .5  km  sweet  spot  position  are  plotted.  For  the 
case  of  Nuri  3,  the  sampling  area  is  biased  to  the  northeast 
of  the  circulation  center  (RLll,  Fig.  7),  and  the  observa¬ 


tional  coverage  within  the  two-degree  radius  circles  is  less 
than  in  Nuri  1  and  2.  This,  combined  with  the  bin  sizes  used 
to  constmct  the  histogram,  leads  to  some  of  the  noisiness 
seen  in  the  Nuri  3  traces.  The  shape  of  the  distributions  in 
the  largest  radius  agrees  fairly  well  with  the  distributions  in 
Fig.  10  of  RLl  1 .  The  most  noticeable  exception  is  that  the 
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Figure  9.  Four-panel  plot  of  relative  frequency  of  absolute  vorticity  for  Nuri  1  (blue),  Nuri  2  (red),  and  Nuri  3  (green)  derived  from  the  10km 
SAMURAI  analyses.  The  relative  frequency  distributions  are  calculated  in  0.5  degree  radius  intervals  from  the  1.5  km  altitude  sweet  spot 
position  out  to  two  degrees.  The  dashed  black  line  represents  the  zero  absolute  vorticity  value. 
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Figure  10.  Two-panel  plot  of  low-level  stretching  tendency  (shaded)  for  Nuri  1  (left)  and  Nuri  2  (right).  The  maximum  stretching  tendency 
between  the  surface  and  1  km  is  displayed.  The  arrows  are  wind  vectors  in  the  co-moving  frame  at  1.5  km  altitude  and  the  thick  black  lines 
are  the  25dBZ  contour  from  the  ELDORA  data  at  5  km.  The  black  box  is  the  four-degree  length  box  centered  on  the  1.5  km  sweet  spot 
position  and  is  for  reference  only.  The  abscissa  is  longitude  and  the  ordinate  is  latitude. 


Nuri  3  distribution  does  not  have  a  large  negative  area  in  our 
plot.  This  is  likely  attributed  to  plotting  the  data  by  radius, 
rather  than  by  the  entire  sampled  area.  Many  of  the  nega¬ 
tive  values  in  RLll’s  distribution  appear  to  be  located  be¬ 
yond  the  two-degree  radius.  This  could  lead  to  a  different 


interpretation  of  the  results  and  illustrates  one  of  the  advan¬ 
tages  in  averaging  over  discrete  radial  intervals.  In  Fig.  9,  the 
distribution  maximum  for  Nuri  1  remains  in  essentially  the 
same  location  at  each  radius;  however,  the  spread  in  the  dis¬ 
tribution  decreases  at  the  0.5  degree  radius.  For  Nuri  1,  this 
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Figure  11.  Two-panel  vertical  cross  section  of  the  frequency  of  pos¬ 
itive  OW  values.  The  abscissa  is  radial  distance  from  the  1 .5  km  de¬ 
rived  sweet  spot  position  in  degrees  and  the  ordinate  is  height  (km). 
The  thick  black  line  represents  the  0.5  contour  (i.e.,  the  location 
where  the  frequency  of  occurrence  of  positive  OW  is  equal  to  the 
frequency  of  occurrence  of  negative  OW). 


suggests  a  somewhat  uniform  voiticity  distribution  through¬ 
out  the  wave  pouch  with  no  strong  maxima.  The  spread  of 
the  Nuri  2  distribution  also  decreases  with  decreasing  radius. 
This  suggests  that,  over  time,  low-level  cyclonic  vorticity  is 
concentrated  near  the  sweet  spot  position.  By  the  Nuri  3  mis¬ 
sion,  the  vortex  is  well  organized,  and  this  is  represented  by 
relatively  large  values  of  low-level  cyclonic  absolute  vortic¬ 
ity  in  the  inner- most  radius. 

Figure  10  shows  the  low-level  stretching  tendency  for 
Nuri  1  and  2  overlaid  with  wind  vectors  in  the  co-moving 
frame  and  the  5  km  25  dBZ  contour.  The  stretching  tendency 
is  calculated  following  RLll,  where  the  divergence  of  the 
advective  flux  of  vorticity  is; 

-V  •  Zi  = -Uft  •  Vh??j  -  (9) 

The  first  term  on  the  LHS  of  Eq.  (9)  is  the  advection  of  abso¬ 
lute  vertical  vorticity  in  the  horizontal,  a  process  that  merely 
redistributes  the  vorticity  in  the  horizontal  plane.  The  sec¬ 
ond  term  is  typically  referred  to  as  the  stretching  term.  Ver¬ 
tical  stretching  (compression)  of  the  column  can  lead  to  a 
localized  increase  (decrease)  in  the  vorticity  magnitude.  Ar¬ 
eas  of  maximum  positive  stretching  tendency  are  observed 
within  the  wave  pouch  near  13.4°  N,  145.8°  E  in  Nuri  1  and 
14.5°  N,  140.1°  E  and  16.8°  N,  140.3°  E  during  Nuri  2,  and 
these  maxima  are  associated  with  cumulus  convection  as  ev¬ 
idenced  by  their  approximate  co-location  with  high  radar  re- 
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Figure  12.  Four-panel  height/di  stance  cross  section  of  each  term 
from  the  circulation  tendency  equation  (Eq.  7)  for  Nuri  1 .  The  ab¬ 
scissa  is  the  length  of  the  box  (degrees)  on  which  the  integration  is 
performed  and  is  centered  on  the  1.5  km  sweet  spot  position.  The 
ordinate  is  height  (km). 


flectivity  at  mid-levels  (Eig.  10).  Infrared  satellite  imagery 
(not  shown)  verifies  the  presence  of  deep  convection  in  these 
regions  at  these  time  periods  (see  Eig.  10, 00:00  UTC  16  Au¬ 
gust,  in  Montgomery  et  al.,  2010a).  Co-location  of  the  max¬ 
ima  in  areas  of  low-level  stretching  tendency  with  deep  con¬ 
vection  suggests  concentration  of  cyclonic  vorticity  through 
vortex  tube  stretching  within  the  wave  pouch.  The  dynam¬ 
ical  processes  observed  here  are  consistent  with  the  mecha¬ 
nisms  necessary  to  develop  a  tropical  cyclone  strength  vortex 
as  outlined  in  DMW09  and  illustrated  by  way  of  explicit  nu¬ 
merical  simulation  in  Montgomery  et  al.  (2010b)  and  Wang 
et  al.  (2010a,  b)  at  cloud-permitting  scales. 

To  assess  the  possible  organization  and  upscale  vor¬ 
ticity  growth  near  the  sweet  spot,  we  will  examine  a 
height/distance  cross  section  of  the  frequency  of  occurrence 
of  OW  values  greater  than  zero  (Eig.  11).  The  thick  black 
line  highlights  the  50  %  contour.  Nuri  1  exhibits  a  layer  of 
greater  than  50  %  positive  OW  from  the  sweet  spot  out  to 
approximately  1.5  degrees  radius.  This  layer  of  positive  OW 
extends  vertically  to  4  km.  Nuri  2  shows  that  while  the  radius 
of  greater  than  50  %  coverage  of  positive  OW  does  not  ex¬ 
tend  much  farther  outward  from  the  sweet  spot  position,  the 
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depth  of  this  positively  dominated  OW  area  now  extends  con¬ 
tiguously  up  to  7.5  km.  At  radii  within  0.6  degrees  from  the 
sweet  spot  position,  the  frequency  of  positive  OW  is  greater 
than  80  % .  These  findings  are  consistent  with  the  building  of 
the  cyclonic  vorticity  monolith  near  the  sweet  spot  through 

(i)  vorticity  intensification  near  the  sweet  spot  caused  by  vor¬ 
tex  stretching  associated  with  increased  convection  between 
Nuri  1  and  2;  (ii)  increasing  longitudinal  alignment  of  the  cir¬ 
culation  in  the  vertical;  and  (hi)  aggregation  of  vorticity  near 
the  sweet  spot  in  the  low-  and  mid-levels  (above  5  km). 

Let  us  now  move  back  up  in  horizontal  scale  in  order 
to  assess  the  spin-up  of  the  pre-Nuri  disturbance  on  the 
mesoscale.  Figure  12  shows  the  analyzed  circulation  ten¬ 
dency  (Eq.  7)  for  Nuri  1  using  the  10  km  SAMURAI  analysis 
around  boxes  of  progressively  increasing  size  in  0.2  degree 
length  increments  centered  on  the  1 .5  km  sweet  spot  posi¬ 
tion.  The  top  row  of  Fig.  12  shows  the  net  tendency  of  the 
absolute  circulation  during  Nuri  1  over  the  analyzed  area  in  a 
distance/height  format.  The  next  three  rows  show  the  individ¬ 
ual  terms  comprising  the  circulation  tendency  equation  (Eq. 
7)  in  the  same  distance/height  format.  The  results  indicate 

(i)  Low-level  (0  to  4  km  altitude)  spin-up  tendency  out  to 
0.8  degrees  from  the  sweet  spot.  In  this  interior  region, 
the  spin-up  is  a  result  of  the  influx  of  absolute  cyclonic 
vorticity  dominating  frictional  spin-down. 

(ii)  Mid-tropospheric  (4-7  km  altitude)  spin-down  tendency 
from  the  sweet  spot  to  approximately  two  degrees  dis¬ 
tance.  Vortex  misalignment  plays  a  role  in  creating  this 
area  of  mid-level  spin-down  tendency,  as  seen  in  the 
negative  convergence  tendency  at  mid-levels. 

(iii)  Barring  a  shallow  region  (<  1  km  depth)  near  the  sur¬ 
face  for  2  and  2.7  degree  boxes  around  the  sweet  spot, 
there  is  a  clear  spin  up  tendency  throughout  the  ob¬ 
served  troposphere  for  boxes  greater  than  2  degree 
lengths. 

(iv)  Maximum  frictional  spin  down  tendency  occurs  in  the 
boundary  layer  beyond  2.25  degree  lengths. 

(v)  An  area  of  low-level  spin-down  tendency  located  within 
the  1-2  degree  boxes  from  the  sweet  spot  (connecting  to 
the  aforementioned  surface-based  spin-down  tendency). 

While  the  surface-based  spin-down  found  for  the  2-2.7  de¬ 
gree  lengths  is  primarily  caused  by  frictional  spin-down  ex¬ 
ceeding  the  influx  of  absolute  vorticity,  the  region  of  spin- 
down  above  the  boundary  layer  between  1-2  degree  box 
lengths  is  a  combination  of  friction,  negative  eddy  fluxes  (not 
shown),  and  negative  contributions  from  the  tilting  term.  In 
summary:  the  foregoing  analysis  suggests  that  there  is  a  com¬ 
bination  of  low-level  spin-up  (where  convergence  of  abso¬ 
lute  vorticity  dominates  frictional  spin-down)  and  spin-down 
within  the  wave  pouch.  However,  on  the  system  scale,  these 
data  indicate  that  the  pre-Nuri  disturbance  during  the  first  day 


of  observations  is  not  spinning  down,  but  broadly  undergoing 
spin -up  in  the  low  levels. 

3.3  System-scale  view  of  Nuri’s  early  spin-up 

A  complementary  way  to  investigate  whether  the  system- 
scale  circulation  is  spinning  up  may  be  obtained  by  exam¬ 
ining  the  time  evolution  of  azimuthally  averaged  absolute 
angular  momentum  surfaces  between  Nuri  1  and  2  (e.g., 
Montgomery  and  Smith,  2014,  their  Fig.  13).  In  this  pa¬ 
per,  Fig.  13a,  b  shows  radius-height  cross  sections  of  az¬ 
imuthally  averaged  tangential  wind  (top,  shaded)  and  abso¬ 
lute  angular  momentum  (top,  contours)  for  the  Nuri  1  (a)  and 
Nuri  2  (b)  flights,  respectively.  Figure  13c  (bottom)  shows 
the  difference  in  azimuthally  averaged  tangential  wind  be¬ 
tween  Nuri  2  and  Nuri  1 .  The  absolute  angular  momentum  is 
defined  by  M  =  n;  -I-  fr^jl,  where  r  is  the  radius  from  the 
sweet  spot  position,  v  is  the  azimuthally  averaged,  storm- 
relative  tangential  velocity,  and  /  is  the  Coriolis  parameter. 
Above  the  frictional  boundary  layer,  M  is  approximately  ma¬ 
terially  conserved.  Therefore,  as  rings  of  fluid  are  drawn  in¬ 
wards,  the  corresponding  tangential  wind  increases. 

In  the  case  of  Nuri  1,  the  maximum  tangential  wind  is 
in  the  lower  levels  and  is  located  near  the  outermost  ranges 
of  the  domain.  For  Nuri  2,  the  maximum  low-level  tangen¬ 
tial  wind  has  moved  inward,  closer  to  the  sweet  spot,  and 
the  magnitude  of  the  tangential  wind  has  increased  in  the 
lower  and  mid- troposphere.  The  increase  of  tangential  wind 
in  the  mid-troposphere  is  consistent  with  an  improved  verti¬ 
cal  alignment  of  the  vortex  as  discussed  above  and  suggests 
a  deepening  wave  pouch. 

The  M  surfaces  show  higher  values  of  M  moving  inwards 
between  Nuri  1  and  2.  As  discussed  above,  inward-moving 
M  surfaces,  in  conjunction  with  its  approximate  material 
conservation  above  the  boundary  layer,  imply  an  amplifica¬ 
tion  of  the  pouch-relative  tangential  wind  field.  These  char¬ 
acteristics  are  consistent  with  a  system  undergoing  develop¬ 
ment  between  Nuri  1  and  2  and  affirm  that,  on  the  system 
scale,  the  pre-Nuri  disturbance  is  broadly  undergoing  spin- 
up  in  the  low  levels  between  Nuri  1  and  2. 


4  Discussion:  relation  with  prior  work  of  RLll  and 
Raymond  et  al.  (2011) 

The  present  kinematical  and  dynamical  analysis  of  Nuri’s 
genesis  observed  during  the  TCS-08  experiment  has  pro¬ 
duced  some  very  similar  results  to  those  found  in  RLll  and 
Raymond  et  al.  (2011).  In  particular,  we  have  found  similar 
distributions  of  low-  and  mid-level  vertical  vorticity;  similar 
areas  of  positive  vortex  stretching  in  the  PBL  and  flow  inte¬ 
rior  that  are  co-located  with  deep  moist  convection;  a  simi¬ 
lar  southeastward  vortex  tilt  between  the  low-  and  mid-level 
mesoscale  circulations;  and  a  similar  vertical  wind  shear  pat¬ 
tern  (not  shown).  The  consistency  of  the  major  structural 
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Figure  13.  Three-panel  plot  of  tangential  wind  (top ,  shaded)  and  ab¬ 
solute  angular  momentum  (top,  contours)  for  Nuri  1  (a)  and  Nuri  2 
(b).  Angular  momentum  surfaces  are  plotted  at  0.2  x  10®  m^  s“^ 
contours  (black)  and  are  highlighted  (white)  at  0.5  x  10®m^s“* 
intervals.  The  difference  in  tangential  wind  speed  between  Nuri  2 
and  Nuri  1  is  plotted  in  (c).  The  abscissa  is  radial  distance  (degrees) 
from  the  1 .5  km  sweet  spot  position.  The  ordinate  is  height  (km). 


features  of  Nuri’s  genesis  using  independent  analysis  tech¬ 
niques  and  software  suggest  that  our  results  are  not  overly 
sensitive  to  the  analysis  methodology. 

The  main  difference  between  these  studies  is  that  RLll 
show  that  the  circulation  tendency  below  approximately 
1.25  km  is  negative  and  thus,  they  conclude  that  the  pre- 
Nuri  circulation  is  decaying  in  the  PEL  during  the  Nuri  1 
mission.  On  the  other  hand,  our  results  indicate  that  the  low- 
level  circulation  on  the  system  scale  is  not  spinning  down  at 
all  radii* .  In  fact,  a  diagnosis  of  the  circulation  tendency  dur¬ 
ing  Nuri  1  at  varying  distances  from  the  sweet  spot  center  (as 
opposed  to  vertical  profiles  calculated  within  a  single  area  as 
employed  by  RLll)  shows  that  the  pre-Nuri  disturbance  is 
already  spinning-up  in  the  lower  levels  (below  4  km  altitude) 
at  various  radii  from  the  sweet  spot  position.  Moreover,  be¬ 


*Because  the  areas  in  which  circulation  tendency  calculated 
herein  and  that  presented  in  RLll  differ,  the  comparison  between 
the  two  results  is  not  direct.  However,  we  have  gone  to  great  lengths 
to  prove  that  both  analysis  schemes  produce  broadly  similar  results 
(see  Sect.  3.2  and  our  response  to  an  Interactive  Comment).  Ad¬ 
ditionally,  one  of  the  initial  results  we  produced  during  the  early 
stages  in  our  mesoscale  analysis  of  Nuri’s  spin-up  sequence  was  an 
attempt  to  exactly  replicate  Fig.  16  from  RLll  (not  shown).  These 
results  are  similar  enough  to  give  us  confidence  that  comparisons 
between  the  circulation  tendency  analyses  are  appropriate.  Thus, 
although  the  circulation  tendency  comparison  may  not  be  a  precise 
“apples-to-apples”  comparison,  our  methodology  provides  new  and 
useful  information  that  is  important  to  understanding  Nuri’s  gene¬ 
sis  process,  information  which  is  omitted  in  RLl  1  and  Raymond  et 
al.  (2011). 


Figure  14.  Vertical  profile  (a)  of  the  mean  virtual  temperature 
difference  over  the  observational  domain  from  our  SAMURAI 
analysis,  black,  and  the  dropwindsonde  analysis  of  Montgomery 
and  Smith  (2012),  red.  Radius-height  cross  section  (b)  of  the  az- 
imuthally  averaged  virtual  temperature  difference. 


tween  Nuri  1  and  Nuri  2  we  show  that  the  low-level  tangen¬ 
tial  wind  increases  at  all  radii  from  the  sweet  spot  in  both  the 
dropwindsonde  and  SAMURAI  analyses.  The  methodology 
of  anchoring  our  analysis  to  the  sweet  spot  position  is  op¬ 
timally  suited  to  observe  and  quantify  the  development  that 
ensues  around  the  translating  sweet  spot,  the  attractor  point 
predicted  by  the  new  cyclogenesis  model  (Montgomery  et 
al.,2012). 

We  now  consider  two  additional  quantities  that  will  allow 
us  to  (i)  further  test  the  thermodynamic  control  hypothesis 
presented  in  Raymond  et  al.  (201 1)  and  RLl  1  and  (ii)  explain 
the  differences  in  the  findings  they  present  compared  to  the 
findings  presented  herein. 

4.1  Mean  virtual  temperature  profiles 

First,  we  examine  thermodynamic  data  from  the  SAMURAI 
analysis  in  order  to  assess  quantitatively  the  thermodynamic 
changes  that  have  been  hypothesized  to  be  a  critical  influence 
on  the  spin-up  process.  Figure  14a  is  a  vertical  profile  of  the 
virtual  temperature  difference  between  Nuri  1  and  2,  aver¬ 
aged  over  the  observational  domain.  Aside  from  the  smooth¬ 
ness  of  the  SAMURAI -derived  profile,  on  the  pouch  scale, 
there  is  little  difference  between  the  profile  derived  from  the 
SAMURAI  analysis  and  that  derived  by  Montgomery  and 
Smith  (2012),  who  used  only  dropwindsonde  data.  These 
data  suggest  only  a  modest  mid-tropospheric  warming  of  1 
K  and  a  very  slight  cooling  (<  0.5  K)  in  the  troposphere  be¬ 
low  2.75  km. 
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A  more  detailed  look  at  the  thermodynamic  structure 
within  the  wave  pouch  may  be  obtained  by  examining  the 
azimuthally  averaged  evolution  of  virtual  temperature  as 
a  function  of  radial  distance  from  the  pouch  center.  This 
methodology  allows  for  an  assessment  of  whether  there  are 
regions  within  the  wave  pouch  that  exhibit  increased  thermal 
stabilization  and  permits  one  to  test  the  statement  by  Ray¬ 
mond  et  al.  (2011;  p.  12)  that  the  “. .  .cold  core  at  low  lev¬ 
els  was  significantly  enhanced  in  Nuri2  compared  to  Nuril”. 
Figure  14b  shows  a  radius-height  cross  section  of  the  change 
in  virtual  temperature  from  Nuri  1  to  2.  The  figure  shows  that 
the  strongest  mid-level  warming  occurs  near  the  center  of 
the  disturbance,  and  the  low-level  cooling  is  extremely  weak 
(less  than  0.5  K),  shallow,  and  confined  largely  to  radii  out¬ 
side  of  1  degree.  In  other  words,  we  do  not  see  evidence  that 
the  cold  core  was  “significantly  enhanced”  between  Nuri  1 
and  2.  Consistent  with  Montgomery  and  Smith  (2012)  and 
Smith  and  Montgomery  (2012),  we  interpret  the  mid-  and 
upper-level  warming  as  simply  the  warm  core  of  the  devel¬ 
oping  tropical  depression  vortex  via  thermal  wind  balance. 
Qualitative  consideration  of  thermal  wind  balance  for  a  cy¬ 
clonic  vortex  in  gradient  wind  balance  indeed  affirms  that 
the  broad  amplification  of  the  low-level  cyclonic  winds  in 
Fig.  13c  is  consistent  with  the  warm-core  thermal  structure 
near  the  center  of  the  pouch  in  Fig.  14b.  These  findings 
are  consistent  also  with  the  numerical  simulations  of  Wang 
(2012)  who  shows  that  the  inner  pouch  region  is  character¬ 
ized  by  mid-level  warming  and  that  weak  cooling  is  found 
only  at  larger  radii  during  the  early  tropical  cyclone  forma¬ 
tion  stage. 

The  foregoing  findings  are  consistent  with  Montgomery 
and  Smith  (2012)  in  the  sense  that  there  is  not  a  strong  indi¬ 
cation  of  the  thermodynamic  stabilization  described  by  Ray¬ 
mond  et  al.  (2011)  on  the  pouch  scale.  Moreover,  on  the  sub¬ 
pouch  scale,  the  data  does  not  suggest  a  significant  enhance¬ 
ment  of  a  cold  core  between  Nuri  1  and  2  at  low  levels  above 
the  boundary  layer^ . 

4.2  Vertical  mass  fiux  profiles 

The  second  quantity  we  examine  is  the  vertical  mass  flux**’. 
Figure  15  shows  the  vertical  mass  flux  profiles  calculated 


^Raymond  et  al.  (2011)  use  saturated  specific  moist  entropy  as  a 
proxy  for  temperature  at  a  given  altitude.  Their  Figs.  2  and  3  show 
a  noticeable  decrease  in  this  quantity  at  the  surface  (~  20  J  kg  * ) 
between  Nuri  1  and  2.  At  the  0.625  km  vertical  level,  however,  their 
analysis  suggests  that  the  decrease  in  saturated-specific  moist  en- 
fi'opy  is  significantly  smaller,  less  than  5  J  kg  K“  * .  Finally,  between 
the  1.25  km  and  5  km  levels,  there  is  no  noticeable  difference  in 
saturated-specific  moist  entropy. 

***For  the  vertical  mass  flux  calculation,  we  reproduced  the  main 
results  of  RLll  using  our  independent  3-D- VAR  analysis  scheme. 
We  integrated  Equation  27  in  RLll  over  the  observational  area  as 
defined  in  RLl  1  at  each  level  to  produce  vertical  mass  flux  profiles 
(Fig.  15). 


Figure  15.  Vertical  mass  flux  profiles  calculated  over  the  obser¬ 
vational  domain  derived  from  the  SAMURAI  analysis  for  Nuri  1 
(blue)  and  Nuri  2  (red). 


over  the  observational  domains.  The  mass  flux  profiles  ap¬ 
pear  broadly  similar  to  those  presented  in  RLl  1 ,  and  confirm 
the  lowering  of  the  maximum  vertical  mass  flux  in  the  ver¬ 
tical  from  Nuri  1  to  2.  Now,  the  lowering  of  the  mass  flux 
profile  is  consistent  with  a  change  from  a  “top-heavy”  latent 
heating  profile  (representative  of  a  mixture  of  deep  convec¬ 
tion  and  stratiform  convection  processes,  e.g.,  Houze  1989, 
2004;  Mapes  and  Houze  1995)  to  a  profile  more  representa¬ 
tive  of  deep  convective  processes.  However,  there  are  subtle 
differences  between  the  RLll  mass  flux  profiles  and  those 
presented  herein  that  invite  careful  examination.  For  Nuri  1, 
the  mass  flux  in  RLll  (their  Fig.  16)  is  approximately  zero 
from  the  surface  to  about  3  km,  and  is  negative  in  some  of 
this  vertical  region.  The  RLl  1  analysis  indicates  also  that  the 
mass  flux  increases  above  4  km,  with  a  maximum  near  10- 
1 1  km.  On  the  other  hand,  our  Nuri  1  analysis  suggests  that 
the  mass  flux  increases  gradually  from  zero  at  the  surface  to 
2  km,  followed  by  a  more  rapid  increase  until  approximately 
5  km.  Between  5-10  km,  there  is  a  broad  region  of  nearly 
constant  mass  flux,  with  an  overall  maximum  near  9  km,  and 
a  secondary  maximum  near  7.5  km.  For  Nuri  2,  our  mass  flux 
profile  displays  an  overall  maximum  near  4  km,  but  decreases 
a  bit  more  sharply  above  the  4  km  maximum  than  the  profile 
presented  in  RLll. 

The  quantitative  differences  highlighted  above  have  im¬ 
portant  physical  ramifications.  Although  both  analyses  have 
vertical  mass  flux  maxima  in  the  upper  troposphere  for 
Nuri  1 ,  our  analysis  indicates  positive  mass  flux  values  in  the 
lower  troposphere  as  well.  Our  analyses  suggest  also  that  the 
influence  of  the  convective  updrafts  exceeds  that  of  the  strat¬ 
iform  downdrafts  in  the  lower  troposphere.  In  contrast,  in  the 
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RLl  1  analysis  the  near  zero  lower  tropospheric  mass  flux  im¬ 
plies  that  on  the  system  scale,  convective  updrafts  are  essen¬ 
tially  balanced  by  downdrafts  produced  in  the  stratiform  rain 
region.  The  positive  low-level  mass  flux  in  our  analysis^* 
is  consistent  with  our  dropwindsonde  only  and  SAMURAI 
analyses,  both  of  which  indicate  a  spin-up  of  the  lower  and 
mid-troposphere  from  Nuri  1  to  2. 

Based  on  the  foregoing  results  and  discussion,  the  results 
presented  herein  suggest  strongly  that  the  observed  proto¬ 
vortex  is  well  into  the  cyclogenesis  process,  independent  of 
the  hypothesized  ‘thermodynamic  control’  as  summarized 
in  the  Introduction.  While  we  certainly  recognize  the  im¬ 
portance  of  thermodynamics  in  the  tropical  cyclogenesis  se¬ 
quence  (Smith  and  Montgomery,  2012;  Montgomery  and 
Smith,  2012),  we  believe  that  the  dynamical  and  thermo¬ 
dynamic  processes  described  in  DMW09,  Montgomery  et 
al.  (2010b),  Wang  et  al.  (2010a,  b)  and  Wang  (2012)  lead  nat¬ 
urally  to  low-level  spin-up  and  vortex  formation  in  the  Kelvin 
cat’s  eye  recirculation  region  of  the  tropical  wave  without 
the  need  for  thermal  stabilization  associated  with  very  small 
temperature  anomalies  (~  0.5  K  or  less).  To  the  extent  that  a 
small  mid-level  stabilization  does  occur  outside  the  core  of 
the  pouch  over  the  two  days  of  intensive  sampling  of  Nuri’s 
genesis,  our  findings  suggest  that  this  stabilization  is  more  in¬ 
cidental  than  essential.  In  other  words,  our  analysis  suggests 
that  the  proposed  thermodynamic  control  presented  by  Ray¬ 
mond  may  not  be  a  necessary  condition  for  low-level  spin-up 
in  the  case  of  Typhoon  Nuri. 


5  Conclusions 

In  this  paper  we  have  continued  our  study  of  the  genesis  of 
Typhoon  Nuri  (2008).  Here  we  have  examined  the  kinemat¬ 
ics,  dynamics,  and  thermodynamics  of  the  genesis  of  Nuri  us¬ 
ing  aircraft  data  from  the  USAF  Cl 30  and  NRL  P3  for  two 
missions  on  15-16  and  16-17  August.  The  mesoscale  analy¬ 
sis  presented  herein  is  complimentary  to  the  large-scale  anal¬ 
ysis  presented  by  Montgomery  et  al.  (2010a)  and  the  thermo¬ 
dynamical  analysis  of  Montgomery  and  Smith  (2012). 

During  the  first  mission,  the  disturbance  was  a  tropical 
easterly  wave  and  prior  to  the  second  mission  the  disturbance 
had  been  upgraded  to  a  TD.  Visual  inspection  of  the  drop¬ 
windsonde  data  suggests  that  the  disturbance-relative  tangen¬ 
tial  wind  speed  increased  in  the  low-  and  mid-levels  between 
Nuri  1  and  2.  Using  dropwindsonde  data,  we  have  shown  that 
the  magnitude  of  the  low-level  (at  and  below  800  hPa)  tan¬ 
gential  winds  did  increase  from  Nuri  1  to  2  at  all  radii  within 


It  is  a  useful  reminder  to  recall  that  the  low-level  spin-up  dis¬ 
cussed  herein  is  caused  by  the  system-scale  convergence  of  vorticity 
substance  induced  by  the  aggregate  of  deep  convective  clouds  in  the 
pouch  and  the  inflow  generated  thereto.  The  area-averaged  vertical 
mass  flux  profile  is  not  the  cause  of  the  spin-up.  Rather,  it  is  simply 
an  indicator  of  the  nature  of  the  precipitation  on  the  pouch  scale. 


the  wave  pouch.  These  findings  indicate  that  there  was  low- 
level  spin-up  occurring  between  Nuri  1  to  2. 

To  compliment  the  dropwindsonde-only  results,  analyses 
from  a  newly  developed  variational  scheme  that  synthesizes 
aircraft  radar  and  dropwindsonde  data  are  presented.  Storm- 
relative  circulation  tendency,  absolute  angular  momentum, 
and  tangential  wind  are  computed  at  various  distances  and 
heights  from  the  sweet  spot  center  throughout  the  observed 
domain  for  each  flight.  The  advantages  of  performing  such 
calculations  at  several  radii  from  the  sweet  spot  position  in¬ 
clude  the  ability  to  examine  the  dynamical  processes  at  all 
observed  distances  within  the  wave  pouch  and  the  ability  to 
directly  compare  data  from  multiple  research  flights  that  may 
have  different  areas  of  data  collection.  Results  from  the  vari¬ 
ational  analysis  are  consistent  with  the  dropwindsonde-only 
analysis  in  that  they  indicate  that  the  low-  and  mid-level  tan¬ 
gential  winds  increased  in  magnitude  and  the  tangential  wind 
maximum  moved  inwards  from  Nuri  1  to  2.  These  results 
are  consistent  with  the  vorticity  aggregation  process.  Addi¬ 
tionally,  by  examining  these  data  in  this  fashion,  we  are  able 
to  observe  and  quantify  the  role  of  vertical  alignment  in  the 
building  of  the  cyclonic  vorticity  monolith. 

By  calculating  the  circulation  tendency  for  Nuri  1 ,  we  have 
identified  regions  indicating  both  low-level  spin-up  and  spin- 
down  within  the  wave  pouch.  Further  examination  indicates 
that  the  low-level  spin-down  results  from  a  combination  of 
frictional  spin-down,  negative  vortex  tilting,  and  eddy  flux 
divergence  out  of  the  integration  area.  There  were  areas  of 
low-level  spin-up  located  closest  to  the  sweet  spot  position. 

From  a  vorticity  perspective,  vorticity  is  observed  to  am¬ 
plify  most  significantly  near  the  sweet  spot  position,  initially 
at  the  lower  levels  and  increasing  in  intensity  with  time.  Ar¬ 
eas  of  positive  vortex  stretching  in  the  lower  troposphere 
associated  with  moist  convection  lead  to  vorticity  intensifi¬ 
cation  in  the  lower  levels  within  the  pouch  and  suggest  the 
presence  of  rotating  deep  convection.  The  regions  of  vortex- 
tube  stretching  increased  in  areal  coverage  and  magnitude 
between  Nuri  1  and  2.  All  of  these  findings  are  consistent 
with  those  of  RLl  1 . 

On  the  system  scale,  inward  movement  of  absolute  angu¬ 
lar  momentum  surfaces  confirms  the  convergence  of  absolute 
angular  momentum,  which  is  one  of  two  system-scale  mech¬ 
anisms  responsible  for  tropical  cyclone  spin-up  (Smith  et  al., 
2009;  Montgomery  and  Smith,  2014).  Radial-height  profiles 
derived  from  the  SAMURAI  analysis  show  that  the  maxi¬ 
mum  increase  in  tangential  winds  is  found  to  have  occurred 
in  the  low  levels  near  the  top  of  the  boundary  layer.  This 
increase  is  consistent  with  the  hypothesis  of  a  progressive 
boundary  layer  control  during  the  spin-up  process  (Mont¬ 
gomery  and  Smith,  2011).  In  short,  the  foregoing  properties 
are  consistent  with  those  of  an  intensifying  proto- vortex  from 
the  beginning  of  the  observation  period. 

The  analyses  and  interpretations  presented  herein  have 
provided  answers  to  the  three  questions  posed  in  the  In¬ 
troduction.  First,  the  magnitude  of  the  low-  and  mid-level 


Atmos.  Chem.  Phys.,  14, 8795-8812, 2014 


www.atmos-chem-phys.net/14/8795/2014/ 


L.  L.  Lussier  III  et  al.:  The  genesis  of  Typhoon  Nuri  as  observed  during  the  TCS-08  field  experiment 


8811 


circulation  increased  from  Nuri  1  to  2.  Second,  for  the  case 
of  Typhoon  Nuri,  our  analysis  suggests  that  any  thermal  sta¬ 
bilization  of  the  atmosphere  occurring  during  Nuri  2  appears 
incidental  rather  than  essential  to  the  storm  development. 
Third,  the  distinct  advantages  of  examining  the  dynamics  of 
tropical  cyclogenesis  in  a  frame  of  reference  anchored  on  the 
translating  sweet  spot  position  were  demonstrated  for  these 
observational  analyses. 

The  findings  of  this  study  are  consistent  in  some  respects 
to  that  of  Raymond  and  colleagues,  but  differ  in  their  sug¬ 
gested  key  result  and  related  scientific  implication  that  the 
pre-Nuri  disturbance  was  spinning  down  in  the  boundary 
layer  on  the  first  day  of  observations.  The  findings  herein 
support  the  marsupial  model  positing  that  the  Kelvin  cat’s 
eye  circulation  of  the  parent  wave-like  disturbance  provides 
a  favorable  environment  for  intrinsic  convective  vorticity 
organization  and  low-level  spin-up  on  the  mesoscale.  The 
findings  further  support  the  applicability  of  the  marsupial 
paradigm  for  easterly  wave-like  disturbances  in  the  western 
North  Pacific  basin. 
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